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Abstract 

We consider models of supersymmetry which can incorporate sizeable mixing between different generations 
of sfermions. While the mixing is constrained by the non-observation of various flavour changing neutral 
current (FCNC) processes, there exist regions of SUSY parameter space where the effects of such mixing 
can be probed at colliders. In this work, we explore this possibility by focussing on the slepton sector. 
The sleptons are produced through cascade decays in direct neutralino-chargino (X 2 X 1 ) pair production at 
the LHC. The final state is characterized by 3 leptons and missing energy. We probe the lepton flavour 
violating (LFV) vertex originating from decay and identify a distinct and unambiguous combination of 
the tri-lepton final state which include a lepton pair with same flavour and same sign (SFSS) in addition to 
a pair with opposite flavour and opposite sign (OFOS). This combination of tri-lepton final state containing 
both OFOS and SFSS pairs, can not only suppress the SM background but can also differentiate the flavour 
violating decays of X 2 from its corresponding flavour conserving decays. We present results for various signal 
benchmark points taking into account background contributions assuming two luminosity options (100 fb _1 
and 1000 fb -1 ) for LHC Run 2 experiment. 
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I. INTRODUCTION 


The experiments dedicated towards the investigation of flavour physics are considered to be one 
of best indirect ways to establish the existence of new physics (NP). They play an important role 
in constraining the viability of various new physics scenarios, thereby complementing the direct 
collider searches. The effects which give rise to large flavour changing neutral currents (FCNC), 
can also be potentially probed at the colliders. For instance, the possibility of observing a flavour 
violating Higgs decay at the Large Hadron Collider (LHC) was discussed in IIH3J. Further, an 
observation of a 2.5 a excess in the H —> channel by CMS [4] in the LHC experiment has 

generated a lot of interest in this sector and has led to a plethora of analysis wm- The leptonic 
sector in the Standard Model (SM) is also interesting owing to the absence of FCNC. This can be 
attributed to the massless nature of neutrinos in the SM. The observation of neutrino oscillations, 
which consequently led to a confirmation of the massive nature of left handed neutrinos, resulting in 
a non-zero decay rate for rare processes like [i —> e'y. The predicted branching ratio (BR) in the SM, 
however is negligibly small (~ 1CU 40 ) due to the tiny neutrino mass and is beyond the sensitivity 
of the current flavour experiments. There exist several extensions of the SM which contribute to 
rare processes such as /i —> via loops, enhancing the BR substantially to 1(T 13 - 10" 15 and 
expected to be within the reach of the indirect flavour probes. Needless to say, an observation of 
such processes is a definitive signal of the presence of physics beyond the SM. Therefore, looking 
for a signal of lepton flavour violation (LFV) directly or indirectly is a challenging avenue to find 
NP. Following this argument, we explore the possibility of observing lepton flavour violation at the 
LHC. 

There are several models in literature which discuss the possibility of flavour violation in the 
leptonic sector. In the current analysis we focus on the supersymmetric extensions of the SM 
which can possess soft masses having significant flavour mixing in the mass basis of fermions. 
This can lead to new contributions to the BR of rare processes. For instance, soft masses with 
flavour mixing can arise in see-saw extensions of SUSY [T9l - [2Tj and also inspired by SUSY GUT 
}22P26j . Alternatively, introduction of flavour symmetries [23128], models with messenger matter 
mixing in gauge mediated supersymmetry breaking (GMSB) ]29H33j . models with R-symmetric 
supersymmetry [34, j35|, supersymmetric theories in the presence of extra-spatial dimensions [ 36L 
{38] etc. also lead to flavourful soft masses. Scenarios in which mass splitting lead to flavour 
violation have been considered in [33 [391 So]. Such extensions in general lead to flavoured soft 
masses and depending on the parameters can lead to observable rates for the flavour violating 
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decays in the squark and leptonic sector. 

Flavour mixing in the sfermion mass matrices can be probed at the collider by the flavour 
violating decay of a sparticle of flavour (say i) into a fermion of flavour j where j 7^ i. Flavour 
violating decays of sleptons were studied in the context of e + e _ linear collider [UH37]- I n Ref. 
08] the authors studied the possibility of observing CP violation from slepton oscillations at the 
LHC and NLC. At the LHC, the sleptons can be produced either through Drell-Yan (DY) process 
or by cascade decays from heavier sparticles. Subsequent flavour violating decays of sleptons 
produced by DY were studied in 09] 50] while those produced by cascade decays were studied in 
[5M3. Probing LFV through the measurement of splitting in the mass eigenstates of sleptons 
was considered in [oSl - lGOl . In this paper we report on our study of flavour violation in the leptonic 
sector by producing sleptons in cascade decays through pair production of neutralino-chargino at 
the future LHC experiments. 

Starting with MSSM, we write the most general structure for the slepton mass matrix. The 
constraints on the model from the non-observation of flavour violating processes can be expressed 
by working in the mass-insertion approximation (MIA) (1201 JLLj in terms of bounds on the flavour 
violating parameter <5^ i / j as defined in Eq |2] [62]. A non-zero 8 % j also opens up the possibility 
of flavour violating decay as far as collider implications of flavoured slepton masses are concerned. 

Our goal is to probe the flavour violating decay in the case of first two generations in the slepton 
sector in SUSY. In this context strong bounds exist on the flavour violating parameter, coming 
primarily from the non-observation of ^ [63]. There exist regions of parameter space where 

these bounds can be relaxed owing to cancellations between different diagrams contributing to this 
process, thereby giving access to probe LFV at the colliders. 

In this letter we explore this possibility to look for LFV decays considering neutralino-chargino 
pair production in proton-proton collisions, which eventually leads to three lepton and missing 
energy final state. The tri-lepton final state is characterized by the presence of two leptons with 
opposite flavour and opposite sign combination (OFOS). The presence of LFV in the tri-lepton 
final state is ensured by demanding a combination of same flavour same sign (SFSS) lepton pair 
along with the OFOS combination. While an imposition of this SFSS criteria along with OFOS 
has a tendency to decrease the signal, it aids in suppressing the backgrounds due to SM and SUSY 
significantly. 

The paper is organized as follows: In Section [TT] we discuss the model set-up introducing the 
various parameters relevant for the analysis in the framework of a simplified model. Relevant 
regions of parameter space consistent with the flavour constraints and conducive to be probed at 
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the colliders are identified in this section. In Section III we explain our choice of OFOS and SFSS 
combination to extract the signal with a detailed description of the simulation. The results of the 
simulation for the background and the representative points for the signal events are presented. 


In Section IV we show regions of the parameter space which can be probed at the LHC Run 2 
experiment in the near future. We conclude in Section [Vl 


II. MODEL PARAMETRIZATION 


In this section we introduce the basic model set-up and related parameters necessary to describe 
LFV. In order to reduce the dependence on many parameters, we consider a simplified SUSY model 
(SMS) approach with only left handed sleptons, wino and a bino while decoupling the rest of the 
spectrum. The /r term is assumed to be ~ 1 TeV to make the neutralino/chargino dominantly 
composed of gauginos with a very small higgsino component. In this case, the mass of the second 
lightest neutralino and Xi , the lightest chargino, are roughly the same as ~ M 2 , the mass of the 
51/(2) gauginos. The lightest neutralino which is assumed to be the lightest supersymmetric 
particle (LSP) has mass ~ Mi, same as the mass of the U( 1) gaugino. 

For the slepton sector we focus on the flavour violation in the left handed sector making the 
right handed sleptons very heavy and set the left-right chiral mixing in the slepton mass matrix 
to be negligible. For simplicity, we assume only two generations. With these assumptions, the left 
handed slepton mass matrix in the basis If = (eF,ft-F) is given as, 


m 2 = 


*iu 


l ii2 

4 2 


(i) 


where F denotes the flavour basis (SUPER CKM) for the sleptons. In this basis the flavour 
violating parameter S 12 is parametrised as [201 [61] . 
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rrij 

«12 = , F ■ (2) 

Naturally, this flavour violating parameter 5\2 is coupled to the rates corresponding to flavour 
violating rare decays in the first and second generation lepton sector. Hence an upper bound on 
this parameter exists due to non-observations of these rare decays like n —> ey [63.J, fi — e conversion 
[6JJ and /j, —> eee [65] • 

In order to obtain the mass eigenvalues of the sleptons, the matrix in Eq[l] can be rotated into 
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a diagonal form by an angle 8 given by, 


sin 28 = 


2 m? T 

^12 


2 2 ’ 
m ^ 2 — m L 1 


where m 2 , are the eigenvalues. It can be related to the flavour violating parameter 5\2 as, 


&12 = 


sin 2 0(m 2 2 — m 2 Li ) 


2m 2 L 


(3) 


(4) 


where uil = " lL ^'" L ~ . The structure of the mass matrix, Eq[l] allows for the possibility of 
flavour oscillations similar to neutrino flavour oscillations. The probability P(e F —» //) of a flavour 
eigenstate e f decaying into a muon is given by [l3j, 

(Am 2 ) 2 


P(e F ^n) = sin- 2 8 — 


> BR(pj —> fj,), 


4m 2 T 2 + (Am 2 ) 2 
sin 2 20 BR(jl -> /r) for T <C Am 2 , 


(5) 


with Am 2 = m 2 a — m 2 . The above expression can be re-expressed in terms of the parameter S \2 
from Eq{4j Thus the branching ratio for the flavour violating decay, X 2 y e ® e L 1 X? can be 
computed as, 

BR(x 2 -t e fj, x?) = Blfv BR(X2 e e) BR(e ->• ey?) + e /r ( 6 ) 


Here the suppression factor due to flavour violation is given by, 

EiFi, = sm 2 2f)= , (7) 

where Am ±2 = mi 2 — m^. 

As mentioned before, bounds on <5 12 and hence Blfv can be obtained by taking into account 
the experimental upper limit on the BR(^j, — > ey) < 5.7 x 10~ 13 [63] . The higher dimensional 
operator contributing to this process is parametrized as [66], 


Cfv = e-^ e a a p ( AlPl + ArPr ) h -E"^, 


( 8 ) 


where the model dependence is captured by the Wilson coefficients Al,r. The branching ratio for 
this process is then given by [66], 

BR(„ e 7 ) = ^5T (l^l 2 + MrI 2 ) ■ (9) 

In our considered model, Ar = 0, as the right handed sleptons are assumed to be very heavy. 


Al on the other hand receives three contributions due to chargino, neutralino and bino mediated 


diagrams and is given as [66] , 


Al 


812 (ay , ( M 2 \ 
m 2 \ 47 t n \ m 2 L ) 


ay 
f 4 tt 



02 
+ 4vr 



( 10 ) 
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where f U)C are loop factors defined in [66] with a non-trivial mass dependence of related sparticles 
and ay^a .2 are the U(l)y and SU(2) gauge couplings. 

The analysis can be simplified again by choosing the following parametrization for the mass M± 
of the (LSP) X], 


Mi = 



( 11 ) 


which is the relation at the electroweak scale due to unification of gaugino masses at the GUT 
scale. For the sleptons we choose, 


M2 > m,L > Mi. 


( 12 ) 


This relation assumes that the intermediate sleptons in X 2 decay are produced on-shell by requiring 
that they are lighter than the mass of X 2 — M 2 . Under these assumptions, we try to find the 
available range of parameters allowed by existing // —» e 7 constraints as will be discussed later. 
Figjljshows the region in the M 2 — plane for which the conditions in EqJTlJand 12 are satisfied 
(green region). It depicts the region of parameter space which is of interest as far as collider 
implications are concerned as discussed in this paper. 

The blue region shows the parameter space for which BR(fi —> e'y) < 5.7 x 10 -13 is satisfied 
for 612 = 0.01 in the left plot and for 5\2 = 0.02 in the right plot. As expected, due to the smaller 
value of 5\2 , the blue region in the left plot has a larger overlap with the green region as compared 
to the plot in the right, thereby admitting smaller slepton masses. The orange region in both the 
plots shows the parameter space for which BR(fi —> ey) < 5.7 x 10 -13 is satisfied for 5\2 = 0.1. 
We find that there is virtually no overlap with the region which is of interest to us from the view 
of collider searches. 

It would be interesting to estimate the suppression factor Blfv corresponding to the allowed 
region in the M 2 — plane for the values of <612 in Figjl] As seen in EqJTJ the parameter Blfv , 
which determines the rate for LFV, is sensitive to the mass-splitting Am = mi 2 — and rriL- 
Blfv increases with 5\2 which can only be accommodated with a larger m^. Thus smaller values 
of 5 12 are not conducive to generate a large Blfv- Blfv is also inversely proportional to the mass 
splitting Am. However, it cannot increase indefinitely as Blfv A 1, leading to a lower bound on 
Am. Figj 2 ] demonstrates the contours of constant Blfv in the Am — mL plane. We find that for 
612 = 0 . 02 , slepton in excess of 250 GeV are required to get Blfv A 0 . 1 , while being consistent 
with the flavour constraints (overlap of blue and green region) in Fig. [I] 
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FIG. 1: Region satisfying Eq[II] and [12] (green), while the orange regions satisfy the /i —> ey constraint 
for S 12 = 0.1. The blue regions are allowed by the upper bound on BR(/r —>• ey) for S 12 = 0.01 (left) and 
5\2 = 0.02 (right). Units of mass are in GeV. 




Am 


FIG. 2: Contours of Blfv for <5 i 2 =0.01 (left) and <h 2 = 0.02 (right). The horizontal blue line is excluded 
by BR(n —> ey) for 612 = 0.01 (left) and 612 = 0.02 (right). The units of mass are in GeV. 


III. SIGNAL AND BACKGROUND SIMULATIONS 


As mentioned in the introduction, we probe the signal of LFV in slepton decay producing it via 
cascade decays of sparticles which are produced in proton-proton collisions at the LHC. Here we 
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focus on X 1 X 2 production which eventually leads to a tri-lepton final state as, 


pp 


xt 


( 13 ) 


l^X°i, 

where i,j denote flavour indices (e,p). The flavour violating vertex causes the decay of a slepton 


(/j), coming from X 2 decay, in Eq. 13, into a lepton of flavour lj with i 7 ^ j. It is clear from the 
above process that the signature of LFV is the presence of 3 leptons of which 2 leptons are with 
opposite flavour and opposite sign (OFOS) in addition to missing energy due to the presence 
of two LSP and neutrino. The leptons with OFOS originate from X 2 decay while the third lepton 
comes from the xf decay. Thus, following this decay scenario, it is possible to have 8 combinations 
of tri-leptons, each having at least one OFOS lepton pair as, 


e + e + p ;e e p + ; p e + p ',p + e p + 

e + e~p + \e~e + p~ \p + e + p~ ;p~e~p + . (14) 

On the other hand, the pair production of X1X2 w dl a ^ so gi ve r i se t° tri-lepton final state with a 
flavour conserving decay of X2 *- e - X2 ~ l l + l~Xv Note that this flavour conserving decay scenario 
also results in 8 combinations of tri-lepton final state given as 

e + p + p~ ;e~p + p~ ]p~e + e~ ;p + e + e~ 

p + p + p~ -,e~e + e~ ;e + e + e _ ; p~ p~p + (15) 

out of which four combinations of OFOS exist as seen in the first line of Eqjl5| It is clearly 
a potential background corresponding to the signal channel in Eq|l4] and expected to have the 
same rate as signal. However, a closer look at these two final states in Eq{l4] and [15] reveals a 
characteristic feature. For example, in the case of signal, out of the 8 combinations of tri-leptons 
with OFOS combinations, notice that four combinations shown in the first line in Eq|14[ also possess 
a pair of leptons with same flavour and same sign (SFSS) which are absent in the background final 
states, shown in Eq|15| The rest of the states with OFOS combination in Eq[l4] are identical to 
the final states given in Eq|15| We exploit this characteristic feature to extract the LFV signal 
events out of all three lepton events including all backgrounds. Thus our signal is composed of 
three leptons having combinations of both OFOS and SFSS together, which is an unambiguous 
and robust signature of LFV in SUSY. Note that while choosing a clean signature of LFV decay 
in SUSY, we pay a price by a factor of half as is clear from Eq|14| However, this specific choice of 




combinations in tri-leptons is very powerful in eliminating much of the dominant SM backgrounds 
arising from WZ and tt following leptonic decays of W/Z and top quarks. 

We now discuss our simulation strategy to estimate the signal rates while suppressing the SM 
and SUSY backgrounds. We performed simulations for both signal and background using PYTHIA8 
m at 14 TeV centre of mass energy and applying the following selections: 

• Jet selection: The jets are reconstructed using FastJet pS j and based on anti kx algorithm 
|69| setting the jet size parameter R = 0.5. The jets passing the cuts on transverse momentum 
p J T > 30 GeV, pseudo-rapidity 177 J | < 3.0, are accepted. 

• Lepton selection: Our signal event is composed of three leptons and are selected accord¬ 
ing to the following requirement on their transverse momenta and the pseudo-rapidity: pi /' 2 ' 3 > 
20,20,10 GeV; It/ 1 ’ 2 ’ 3 ! < 2.5, where the leptons are px ordered with p 7 1 being the hardest one. 
In addition, the leptons are also required to be isolated i.e. free from nearby hadronic activities. 
It is ensured by requiring the total accompanying transverse energy, which is the scalar sum of 
transverse momenta of jets within a cone of size A R(l,j) < 0.3 around the lepton, is less than 10% 
of the transverse momentum of the corresponding lepton. 

• Missing transverse momentum: We compute the missing transverse momentum by carrying 
out a vector sum over the momenta of all visible particles and then reverse its sign. Since p/x is 
hard in signal events, so we apply a cut p/x > 100 GeV. 

• Z mass veto: We require that in three lepton events, the invariant mass of two leptons with 
opposite sign and same flavour should not lie in the mass window mu = Mz ± 20 GeV. It helps to 
get rid of significant amount of WZ background. 

• b like jet selection: The b jets are identified through jet-quark matching i.e. those jets which 
lie with in A R(b,j) < 0.3 are assumed to be b like jets. 

• OFOS: Our signal event is characterised by the requirement that it has at least one lepton pair 
with opposite flavour and opposite sign. 

• SFSS: We require the presence of SFSS combination along with OFOS combination in three 
lepton final state, which is the characteristic of our signal. As stated before, this criteria is very 
effective in isolating the background due to the same SUSY process but for their subsequent flavour 
conserving decays, in particular for x® decay. 

We perform our analysis by choosing various representative points in the SUSY parameter space. 
The spectrum is generated using SUSPECT ED] and the decays of the sparticles are computed using 
SUSYHIT m- Table [I] presents the six representative points (A-F) for which we discuss the details 
of our simulation. From A to F, the spectrum is characterized by increasing masses of gauginos, 


9 


with the slepton mass tul lying midway between the two, rriL = (Mi + M 2 ) /2. 


Spectrum Characteristics 

A 

B 

C 

D 

E 

F 

Xi/Xx 

210 

314 

417 

518 

619 

718 

Xi 

95.8 

144 

193 

241 

290 

339 

m L 

156 

229 

303 

377 

452 

526 

BR{xl e L e) 

0.13 

0.15 

0.16 

0.16 

0.16 

0.16 

BR(x 2 -> AlM) 

0.13 

0.15 

0.16 

0.16 

0.16 

0.16 


TABLE I: Representative choices of SUSY parameter space. All masses are in GeV. 


In Table |TT] and III we present the effects of selection of cuts in simulation for both the signal 
and backgrounds respectively. In addition to the SM backgrounds which are mainly due to tt 
and WZ, we also simulate the background taking into account the contributions due to flavour 
conserving decay of x® f° r each of the representative points in Table |l} There are other sub-dominant 
backgrounds like tbW, ZZ if one lepton is missed or WW, if jets fake as leptons. However these 
backgrounds are expected to be very small and not considered here. We present results for signals 
corresponding to those representative parameter space as shown in Table [TT] . In this table, the first 
column shows the sequence of cuts applied in the simulation, while the second column onwards 
event yields for the signal are shown. Table |TTT] presents the same for the backgrounds due to SUSY 
in the second column and the SM in the third column. Notice that lepton isolation requirement 
and a cut on pfp has considerable impact in reducing tt and WZ background. As noted earlier, 
we find the SFSS criteria to be very effective in isolating the SUSY background due to flavour 
conserving decay of X 2 l° r a ll the representative points in Table [I] . Finally, it is possible to have 
large number of tri-lepton events in background processes, but imposition of specific choices like 
OFOS and SFSS along with large missing energy cut help in isolating it to a great extent as shown 


in Table III . In spite of this suppression of background events, the signal yields are far below than 
the total background contribution owing to it’s huge production cross sections as shown in Table 
III Therefore, in order to improve signal sensitivity further, we impose additional requirements by 
looking into the other characteristics of signal events. For example, signal events are free from any 
kind of hadronic activities at the parton level i.e. no hard jets are expected in the signal final state, 
whereas in background process, in particular events from tt are accompanied with large number of 
jets. We exploit this fact to increase signal sensitivity by adding following criteria. 

Case a: Jet veto 


In this case we reject events if it contain any hard jets. In Table III we see that while the jet veto 
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Signal( X ° X ±) 

M 2 => 

200 

300 

400 

500 

600 

700 

No. of events generated 

10000 

10000 10000 

10000 

10000 

10000 

>20 iV t >10, M <2.5 

1371 

1752 

2014 

2218 

2225 

2342 

Lepton isolation cut 

1330 

1669 

1883 

2055 

2036 

2112 

pfr > 100 

474 

959 

1326 

1600 

1683 

1860 

OFOS 

470 

952 

1319 

1581 

1659 

1828 

Z mass veto 

423 

849 

1218 

1485 

1574 

1752 

SFSS 

223 

462 

640 

783 

804 

892 

Case a: jet veto 

91 

205 

288 

337 

346 

380 

Case b: 6-like jet veto 

221 

458 

635 

777 

798 

884 

Case c: rij < 1 and 6-like veto 

161 

375 

479 

604 

617 

687 


TABLE II: Event summary for signal after all selections. All energy units are in GeV. 



SUSY( X ° X ±) 

SM 


A 

B 

C 

D 

E 

F 

tt 

WZ 

m 2 =► 

200 

300 

400 

500 

600 

700 

- 

- 

Cross section (/6) at 14 TeV 

1.65 x 10 3 370.5 

118.8 

45.6 

20.5 

9.57 

9.3 x 10 5 4.47 x 10 4 

No. of events generated 

10000 

10000 

10000 10000 

10000 

10000 

10 7 

3 x 10 6 

p^' 2 > 20,> 10, M < 2.5 

1299 

1779 

2015 

2195 

2245 

2361 

164895 

23960 

Lepton isolation cut 

1251 

1672 

1874 

2044 

2051 

2131 

70233 

22366 

pfr > 100 

454 

967 

1311 

1624 

1722 

1872 

19241 

1669 

OFOS 

209 

482 

656 

820 

855 

918 

14012 

858 

Z mass veto 

126 

346 

547 

728 

768 

853 

12395 

122 

SFSS 

4 

6 

11 

14 

15 

25 

4598 

22 

Case a: jet veto 

< 1 

1 

1 

5 

4 

4 

29 

< 1 

Case b: 6-like jet veto 

4 

5 

10 

14 

13 

23 

131 

13 

Case c: rij < 1 and 6-like veto 

1 

3 

7 

9 

9 

19 

48 

5 


TABLE III: Event summary for SUSY and SM background. All energy units are in GeV. 

criteria reduces the tt and WZ background significantly, but it also substantially damage the signal 
by a factor of 2 or 3 as shown in Table [II]. In signal process, jets arise mainly from the hadronic 
radiation in initial and final states and it is true for all the representative signal points. The reason 
can be attributed to enhancement of hadronic activities at higher energies. Nevertheless the jet 
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veto seems to be useful to improve signal to background ratio. However we consider two more 
alternatives with a goal to increase signal sensitivity further: 

Case b: 6-like jet veto 


Here we eliminate events if there be at least one 6 like jet. As can be seen Table III, 6 jet veto is 
more efficient than the jet veto condition, as the tt background is suppressed by a few orders of 
magnitude without costing the signal too much. 

Case c: Apply 6-like jet veto and number of jets rij < 1 
Here we apply the 6-like jet veto condition along with the presence of maximum one jet. As seen 


in Table III , it is very helpful in reducing the tt background significantly but it does not affect the 
signal as much as the simple jet veto condition (rij = 0) does alone. 

Note that we have identified 6-like jets by a naive jet-quark matching which is an overestinration 
from the realistic b-jet tagging[72] which is out of scope of the present analysis. However, for the 
sake of illustration, we present these results with b-jet veto, (case (b) and (c)), to demonstrate 
that this criteria might be very useful in suppressing backgrounds, which requires more detector 
based simulation. In view of this, we focus only on the results obtained by using jet veto, case(a) 
for further discussion. 

We also present the dilepton (e/r) invariant mass distributions for the spectrum A (left) and 
F (right) in Figj3] normalizing it to unity. It is subject to all primary selection cuts on leptons 
and jets, including the OFOS and SFSS combination. The m eil distribution is expected to have a 
sharp edge on higher side, which can be derived analytically from kinematical consideration. The 
position of this edge of m e/1 is given as [53], [73], 


C m?D 2 = ml o 1 - 



1 - 



(16) 


The appearance of an edge in the m efl distribution is a clear indication of LFV vertex in the 
X 2 decay. However, this m e/J , distribution is affected by a combinatorial problem. For each tri¬ 
lepton event, two OFOS pairs can be constructed: a) both leptons coming from X 2 decay and b) 
“imposter” pair with one lepton from X 2 an d the other from xf- I n Figj3]the red (dotted) curve 
represents the dilepton invariant mass distribution of the leptons tracked to the X 2 vertex while 
blue (solid) curve corresponds to dilepton without any prior information about their origin. It (red 
dotted line) exhibits a very distinct edge as the identity of the lepton pair originating for X 2 is 
known a-priori. The (solid) blue line is more realistic as it includes both the correct OFOS and 
SFSS pair as well as the contamination due to the “impostor” pair which is responsible for a tail 
beyond the edge. As a result it exhibits a more diffused behaviour near the position of the edge. 
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FIG. 3: OFOS dilepton invariant mass distribution for spectrum A (left) and spectrum F (right). The events 
are selected at the SFSS level. 


Properties 

Signal (S) 

Background (B) 

A B C D E F 

it WZ 

Cross section (fb) at 14 TeV 

1.65 x 10 3 370.5 118.8 45.6 20.5 9.57 

9.3 x 10 5 4.47 x 10 4 

Normalized cross sections 

Case a: jet veto 

Case b: 6-like veto 

Case c: rij < 1 and 6-like veto 

15.01 7.59 3.41 1.51 0.67 0.37 

36.4 16.9 7.54 3.54 1.63 0.85 

26.5 13.9 5.7 2.75 1.26 0.66 

2.69 < 1 

12.1 0.19 

4.4 0.07 

^=(@100) fb- 1 

Case a: jet veto 

Case b: 6-like veto 

Case c: rij < 1 and 6-like veto 

91.43 45.93 20.78 9.32 4.31 2.24 

100.99 47.87 21.34 10.04 4.64 2.43 

122.4 64.4 26.4 12.8 5.92 3.12 

- 


TABLE IV: Normalized cross-section (fb) and S/y/B for signal and background subject to three selection 
conditions 

However, we can roughly estimate the position of the edge using the blue (solid) line as ~ 120 GeV 
for the left panel and ~ 375 GeV for the right panel. We find that these values are in fairly good 
agreement with the corresponding numbers used in our simulation. It may be noted here that such 
distributions with a sharp edge are the characteristic feature of these type of decays which can also 
be exploited to suppress backgrounds [53] in order to increase signal to background ratio. 
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IV. RESULTS AND DISCUSSIONS 


Table [TV] gives the normalized signal and background cross-sections due to all selection cuts. 
These are obtained by multiplying the production cross section given in the first row by acceptance 
efficiencies. The production cross section are estimated by multiplying the leading order (LO) cross 
section obtained from PYTHIA8 with the corresponding k factors 1 . Corresponding to these signal 
and background cross sections, we also present the signal significance by computing S/y /B for 


integrated luminosity 100 flu 1 as shown in the bottom of Table IV Although case(b) corresponding 
to b -like jet veto results in the largest cross section for all signal parameter space, signal significance 
does not improve due to comparatively less suppression of SM backgrounds. With the increase of 
gaugino masses acceptance efficiencies goes up as final state particles become comparatively harder, 
but S/y/~B is depleted due to drop in X 2 X 1 pair production cross-section. While estimating signal 
rates and significance, we assume a maximal flavour violation i.e. Blfv = 1- Obviously, a further 
suppression is expected by a factor Blfv which depletes the BR of x/ (see Eqj6|. For a given 
£ 12 , Blfv is a function of the slepton mass as well as the mass splitting A m as shown in Fig. [2j 
For instance S/y/B may suffer by an order of magnitude for Blfv = 0.1. While the lower end of 
the spectrum can lead to a larger S/y/B, the corresponding Blfv decreases as we move further 
towards the IR part of the slepton spectrum. This can be attributed to stronger bounds on 5 12 
for lower slepton masses. Though the lower mass is not yet ruled out, it is more economical to 
consider relatively heavier slepton masses as the bounds from current and future experiments will 
be relatively weaker. 

In FigjiJ we illustrate this mass sensitivity by presenting S/y/B obtained using jet veto condition 
case (a). Notice that for a given xf and X 2 masses, signal is not very sensitive to slepton mass 
as long as it is produced on-shell from X 2 decay and M x 0 — is sufficiently high. The regions 
in the M 2 — niL plane correspond to different values of S/y/B computed for £ = 100 fb _1 and by 
assuming Blfv = 1- The sleptons and gaugino masses follow the parametrisation in Eqjll| and 
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It is superimposed on the region satisfying BR(/j, —>• ey) < 5.7 x 10~ 13 for 612 = 0.01 (left) 
and 8\2 = 0.02 (right). As seen from Table IV and Fig. |4j the signal significance is better for lower 
masses due to larger X 2 X 1 pair production cross section. However, it suffers by smaller values of 
Blfv corresponding to those slepton masses as shown in Fqj7]and Figj2] 

Fig{5] shows the sensitivity reach of Blfv hi the M2 — niL plane using the parametrisation 


1 The appropriate k factors for tt and WZ processes is 1.6 m and 1.7 75] respectively while for the signal it is 

1.5 [7B|. 
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FIG. 4: Variation of S/yfB (using jet veto, case (a) for different regions with two choices of S 12 = 0.01 
(left) and S 12 = 0.02(right). The regions light blue are allowed by BR(/n —>• ey) constraint. Here we assume 
Blfv = I- Masses are in GeV. 


in EqJTl] and 12 The numbers in boxes for different coloured regions give the minimum values 
of Blfv which can be probed, while requiring a 5 a discovery corresponding to those values of 
M 2 and rriL and are presented for two different options of luminosities: £ = 100 fb -1 (left) and 
£ = 1000 fb -1 (right). As the constraints from indirect flavour measurements get tighter, larger 
Blfv can be attained with heavier slepton masses, while respecting bounds from the rare decays 
as shown in Fig. [I] and [ 2 J For example, for lower masses X 2 ~ X\ ~ 250 GeV and ~ 200 

GeV, the LFV parameter Blfv ~ 0-05 or more can be probed at 5cr level of signal sensitivity for 
£ = 100 fb . As expected, the minimum Blfv required for a 5er sensitivity goes up, thereby 
reducing the sensitivity of Blfv measurement with the increase of gaugino and slepton masses and 
this can be attributed to the drop in cross-sections. The left plot in Figj5]is terminated at the point 
corresponding to a requirement of Blfv= 1 f° r a 5u discovery. As a result, the representative points 
E and F corresponding to heavier slepton masses are beyond the sensitivity of LHC at £ = 100 fb 
as they require Blfv > 1 to achieve a 5 a discovery. However, flavour violating decays with heavier 
slepton masses as high as 650 GeV can be probed with an integrated luminosity of £ = 1000 fb 
as shown in the right plot of Fig [5} 
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FIG. 5: Minimum value (in small box) of Blfv for a S/yfB = 5 discovery for £ = 100 fb ~ 1 (left) and 
£ = 1000 fb -1 (right). The S/y/~B is computed using jet veto condition. The filled triangles correspond 
to the representative points A-F from left to right. The plot is truncated at the point where Blfv > 1 is 
required to get a 5 a sensitivity of signal for that particular luminosity. Masses are in GeV. 


V. CONCLUSION 


The observation of flavour violating rare decays would be one of the best indicators of the exis¬ 
tence of physics beyond the SM. Measurements of such decays play an important role in constraining 
several new physics models and hence has received a lot of attention recently. We attempt to ex¬ 
plore the flavour violation in the lepton sector in the context of well motivated models of flavourful 
supersymmetry. We follow an approach based on a simplified model with only the left handed 
sleptons along with the neutralinos which are gaugino dominated. We consider pair production of 
xlxt and their subsequent leptonic decays which includes the LFV decays of x®- The final state is 
composed of three leptons and accompanied by large missing energy. In addition to the presence of 
a lepton pair with OFOS, we observed that certain tri-lepton combinations are also characterized 
by a lepton pair with SFSS -which is a unique and robust signature of LFV in SUSY. 

The discovery potential of observing this LFV signal is primarily dependent on the masses of 
sleptons and gauginos. These masses are however constrained by non-observation of FCNC decays 
such as /r —> ey and they get stronger as the flavour violating parameter 5\2 becomes larger. We 
have identified the allowed range of slepton and gaugino masses relevant for our study. In addition 
variation of LFV parameter Blfv with masses of slepton and mass difference between lepton mass 
eigenstates (Am) are also presented. 
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Estimating the various background contributions, we predict the signal sensitivity for a few 
representative choices of SUSY parameters. The combination of three leptons with OFOS and 
SFSS is found to be very useful to achieve a reasonable sensitivity. It is found that for gaugino 
masses ~ 250 GeV and slepton masses ~200 GeV, the LFV parameter Blfv as low as 0.05 can 
be probed with 100 fb -1 integrated luminosity. For heavier masses ~ 600 — 700 GeV, because of 
reduced X 2 X 1 pair production cross section, the measurement of LFV parameter Blfv requires 
higher luminosity ^1000 fb -1 . Our study clearly establishes the prospects of finding LFV signal 
in this SUSY channel at the LHC Run 2 experiment with high luminosity options. 
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